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We study transport properties of the helical edge channels of a quantum spin Hall (QSH) insula¬ 
tor, in the presence of electron-electron interactions and weak, local Rashba spin-orbit coupling. The 
combination of the two allows for inelastic backscattering that does not break time-reversal sym¬ 
metry (TRS), resulting in interaction-dependent power law corrections to the conductance. Here, 
we use a non-equilibrium Keldysh formalism to describe the situation of a long, one-dimensional 
edge channel coupled to external reservoirs, where the applied bias is the leading energy scale. 
By calculating explicitly the corrections to the conductance up to fourth order of the impurity 
strength, we analyse correlated single- and two-particle backscattering processes on a microscopic 
level. Interestingly, we show that the modeling of the leads together with the breaking of Galilean 
invariance has important consequences on the transport properties. Such breaking occurs, because 
the Galilean invariance of the bulk spectrum transforms into an emergent Lorentz invariance of the 
edge spectrum. With this broken Galilean invariance at the QSH edge, we find a contribution to 
single particle backscattering with a very low power scaling, while in the presence of Galilean in¬ 
variance the leading contribution would be due to correlated two-particle backscattering only. This 
difference is further reflected in different values of the Fano factor of the shot noise, an experimen¬ 
tally observable quantity. The described behaviour is specific to the Rashba scatterer, and does not 
occur in the case of backscattering off a time-reversal breaking, magnetic impurity. 


I. INTRODUCTION 

Quantum spin Hall systems 1-3 are prominent exam¬ 
ples of a topologically non-trivial phase of matter. Even 
though the bulk system is gapped, gapless edge states 
with a linear energy dispersion are formed at the bound¬ 
aries of the (two-dimensional) QSH sample, in con¬ 
tact with vacuum or a normal insulator. These one¬ 
dimensional edge channels are responsible for peculiar 
transport properties and contribute to the conductance 
with Go = e 2 /h per channel. Such quantized values of 
the conductance have been experimentally observed e.g. 
in HgTe/CdTe 4 or InAs/GaSb quantum wells 5 . 

The ballistic edge transport in the QSH system is pro¬ 
tected by two mechanisms - the helical character of the 
transport channels and time-reversal symmetry (TRS). 
First, helicity is an essential feature of the QSH edge 
states. The electrons direction of motion is strongly cou¬ 
pled to the spin direction, meaning that electrons prop¬ 
agating in opposite directions are expected to have op¬ 
posite spins. In inversion-preserving geometries, we can 
describe the edge electrons by the two orthogonal eigen¬ 
states of the s z spin operator. With spin-orbit coupling 
(SOC) present, the spin s z component is not preserved 
any more, however, Kramers theorem ensures the exis¬ 
tence of a pair of spin-orthogonal states with opposite 
momentum, such that the helical character remains 6 . By 
this virtue, the helical edge transport is insensitive to dis¬ 
order that does not flip the electrons spin, since there is 
simply no adequate propagation channel for the backscat- 
tered electron. In the presence of spin-flipping impurities, 
time-reversal-symmetry still prevents elastic backscatter¬ 
ing between counter-propagating channels 1 ’'’ 8 . 
Importantly however, these mechanisms do not pro¬ 


tect the helical edge against inelastic backscattering. 
Coulomb interactions or phonons will lead to single- 
(SPB) or multiple-particle backscattering 6,9-12 . Espe¬ 
cially, the two-particle backscattering (TPB) processes 
can be a relevant perturbation for strong interaction 
strengths, and open up a gap in the spectrum, while SPB 
contributions were previously shown to be vanishing 11 . 
Such findings have been achieved essentially with the 
help of renormalization group (RG) approaches, using 
finite temperatures at equilibrium 9,10,13 . In Ref. 12, a 
non-equilibrium framework was used to study the cor¬ 
rection to the conductance at finite voltage bias due to 
generic, inelastic backscattering. 

In this work, we study both single and two-particle 
backscattering off a single, local Rashba scatterer, start¬ 
ing from a microscopic picture. To account for a more 
realistic situation, we apply a voltage bias on two non¬ 
interacting leads connected to an interacting helical Lut- 
tinger liquid and use a non-equilibrium Keldysh formal¬ 
ism to calculate the expected backscattering current up 
to fourth order of the impurity strength. Somewhat sur¬ 
prisingly, the modeling of the external reservoirs as well 
as the rigorous treatment of elementary invariances, such 
as the Pauli principle and Galilean invariance, have qual¬ 
itative effects on the resulting current. In particular, we 
discover a not yet reported single-particle contribution, 
that may arise in the absence of Galilean invariance and 
the aforementioned design of the leads. The correspond¬ 
ing correction to the conductance is found to scale with 
a remarkably low power of 5G ~ V 2 in the limit of weak 
interactions. On the other hand, if Galilean invariance 
were present, the lowest order contribution would stem 
from correlated two-particle processes. Besides the dif¬ 
ferent power law scalings of the correction to the con- 
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ductance, the distinct character of the two situations is 
also reflected in the Fano factor, as we verify calculat¬ 
ing the shot noise. In the case of broken Galilean in¬ 
variance, we predict a Fano factor of one, while in the 
presence of Galilean symmetry, correlated two-particle 
backscattering processes become manifest in a Fano fac¬ 
tor of two. We therefore present an experimentally ac¬ 
cessible quantity, the shot noise, as a direct evidence of 
broken Galilean invariance in this setup. For a better un¬ 
derstanding of the problem, we furthermore work out the 
differences with the case of a TRS-breaking impurity. 
The article is organized as follows: In Sec. II, we give 
a detailed description of the model, deriving the effec¬ 
tive operators in the presence of the external bias. While 
electron-electron interactions are treated exactly in the 
framework of bosonization, we perform a perturbative 
expansion in the Rashba impurity strength. In Secs. Ill 
and IV, the averaged backscattering current, as well as 
the shot noise, are calculated perturbatively in second 
and fourth order of the impurity strength, respectively. 
In order to verify our findings in Sec. III.A, an alter¬ 
native calculation in fermionic language is presented in 
Sec. III.B. Finally, in Sec. V, we summarize our results. 
Additional details are given in the Appendices. 


II. MODEL 


A. Bosonization 


The system under consideration can be described by a 
Hamiltonian consisting of the four terms H = Hq + Hr + 
H v + H int , with 


H 0 = dx Y tfr*(a) (- irv F d x - p r ) ^(z), 

J r =± 


h r = 


J dx a(x) [(d,^) (d x <F_)] (x) +H.c., 


H V = — 


dx (cc)\R_ 4 _ (a;) + 

Hint =92 J dx ^(®)*l(s)tf_(s)tf+(a;). 


and H , ±(a;) are fermionic creation and annihila¬ 
tion operators for a right (+) or left (—) moving particle 
and vf the Fermi velocity. /i + and /i_ are the chemi¬ 
cal potentials of the right and left moving particles, re¬ 
spectively. We set H = 1 in this article unless explicitly 
stated. H 0 describes the free Hamiltonian with a strictly 
linear dispersion relation 14 and H lri t embodies electron- 
electron interactions, where we take into account only 
processes of the type p 2 coupling right and left movers. 
It is a particular feature of the helical liquid that in the 
case of contact interactions the terms g\ and <74 of the 
usual g-ology are absent, as we explain in the Appendix 
A. The system is interspersed by perturbations that we 
choose to be of the form of a weak, Rashba-like impurity 


H r coupling right and left movers 9,10,13 . a(x) may be 
any function that models the presence of Rashba impu¬ 
rities, however, in this article, we will restrict ourselves 
to a Dirac-like 5-function for simplicity. This impurity 
Hamiltonian Hr constitutes, from our perspective, the 
simplest model that can couple the two edge states and 
preserve TRS at the same time. Note that it has also 
been derived on the basis of rather general assumptions 
in the context of the so-called generic helical liquid 6,15 . 
In the following, we will use the technique of bosonization 
to treat interactions exactly, while including the impurity 
strength perturbatively. We make use of the bosonization 
identity 14,16 - 18 

(x) = F ± —^e ±i ‘^ N±x e- i{± ^ x) - 9 ^ ) , 

V27ra 


where F± is the Klein factor for a right/left moving parti¬ 
cle, N± the particle number operator counting right/left 
moving particles, a a short-distance cutoff, and <fi and 
6 two bosonic fields obeying the commutation relation 18 
[<p(x),d X '6(x')] — in(6(x — x') — j). In the limit of a 
large wire length, L —>• 00, the expressions will simplify, 
however, this limit has to be taken with care. The fields 
cf> and 9 are related to the density and current operators 
respectively, through 

■ p+{x): + : p_( x) := y{N + + N_) — -d x <j>{x) , 

L 7r 

: p+(x) : - : P-{x) := y(N+ - AL) + -d x 9(x) . 

L 7r 

Here, we adopt the definition of the bosonic fields of 
Ref. 14, except for that the zero-modes N + and AT are 
written explicitly and are not included in the fields <f> 
and 9 anymore. We then specify the chemical poten¬ 
tial introducing the voltage bias V by p± = ±eV/2. 
The Hamiltonian, in its bosonized form, is now given 
by H = H ll + H v + Hr , with 10 


H LL = ^(Nl + N*) 


V 

27r . 


dx 


K:(d x 6(x)) 2 : + — :(d x (f)(x)) 2 : 


eV 

H V = -(N + -N_), 


Hr = iF\_F- / dx 


a(x) / 2 ™ 
ira \ L 


K 


(1) 

( 2 ) 


x : (d x B{x) + ^(N+- iV_)) e MH^)-^(N + +N_)) : 
+ H.c. (3) 


Operators between colons :(...): are normal-ordered. 
The free Luttinger liquid Hamiltonian Hll comprises 
both the linear dispersion and the Coulomb interac¬ 
tions, embodied in the new sound velocity 14 v = 
vfO- — ( 2 ^ ) 2 ) 1 / 2 and the dimensionless parameter K = 

( i+ 32 /( 2 g«y) ) 1 ^ 2 • non-interacting case <72 —> 0 cor¬ 
responds to A' —► 1, while 0 < K < 1 means finite, 
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repulsive interactions. From Eq. (1), one can derive the 
bosonic form of the continuity equation 

d t 4> - vKd x 9 = 0 . (4) 

The fact that vK = Vf{1 — 92 /(27rvp)) appears in this 
equation, instead of vp, is a hallmark of ID interacting 
Dirac fermions. Density operators actually do not 
commute, leading to a renormalization of the current 
operator by interactions. 

With this setup in mind, we now follow the lines of 
Ref. 19, with the zero modes unaffected by electron- 
electron interactions. We emphasize that this point 
will have important consequences on the transport 
properties, as we show in the following. The idea behind 
this model is that the physical contacts attached to the 
quantum spin Hall edge are usually higher-dimensional, 
metallic objects of a greater extent than the edge channel 
itself. Since screening can work more efficiently here, we 
expect the leads to be barely interacting 20-22 . The full 
system is a composition of two external leads and the 
edge channel in between (see Fig. 1). Thus, in the limit 
of large lead sizes, the zero modes of the full system 
will be dominated by the zero modes of the contacts 19 . 
Because of that, the influence of the non-interacting 
leads can be incorporated in the decoupled zero-modes. 
In a second step, we then identify the length of the edge 
channel itself, L, as much greater than the remaining 
lengthscales of the system, implying the hierarchy 
L Vp/eV a. Taking into account that position 
integrals can potentially compensate factors of 1 /L, we 
can apply the limit L —» 00 to eventually drop the zero 
modes (and Klein factors). 



Figure 1. (Color online) Scheme of the interacting QSH edge 
channel, attached to two long, non-interacting leads. The 
external bias modifies the number of right and left moving 
electrons in the contacts. Transport along the edge channel is 
interrupted by a local Rashba SOC scatterer (dot). The inset 
shows the linear energy spectrum of the edge state electrons. 

ing current or the current-current correlation (noise). If 
a non-zero bias voltage is applied, this will modify the 
expectation value of any observable O , for instance the 
current, to 19 

(O) = 4 Tr (e- 0 H ' v e iHt Oe- lHt ). (6) 

Zj 


To explore the underlying mechanisms, it is interesting 
to compare our findings to the well known problem of a 
regular impurity 23 . Since at the quantum spin Hall edge 
backscattering is only possible with a simultaneous spin 
flip, the Rashba impurity Hp in this case is replaced by 
an impurity of the form 


H m = 


dx m(x) ^\E , + (:t)\E , _(:e) +H.c.^ 


K 


X:e 2iWx) = * ? (iV + +iV_)) ; + Hx _ ) 


( 5 ) 


with impurity strength m(x). Such an impurity can 
cause backscattering of a right-moving particle into a 
left-moving particle, however, it does not preserve time- 
reversal symmetry, which is why we refer to it as a mag¬ 
netic impurity. Both perturbations will hinder electronic 
transport in a characteristic form. 


B. Non-equilibrium transport 

Transport signatures of the system are provided quite 
generally by quantities such as the average backscatter- 


Here, we have defined Z = Tr(e -/lff v ), as well as 
H = Hll + Hr and H'y = Hll~ Hy. As we describe in 
Appendix B, the external bias can be implemented most 
conveniently by a unitary transformation U on the free 
system, resulting in a shift of the bosonic fields. The ex¬ 
pectation value of the observable O due to backscattering 
is then given by 19 

(O) = Tr(e~^SHt) (o) (t)S(i)), (7) 

All v 'i 


with 


S(t) = T exp 


—1 


df (H LL + [H R y ? (o) 


The operator labels denote O = WOU and Oj(t) = 
e iH v tQ e -iH v t j n {.jjg interaction picture. U is the unitary 
operator specified in Eq. (Bl) and Zpp = Tr(e -/3ffiI '). 
The effect of external bias is thus fully absorbed in the 
two effective operators (£) an d (t). 

To calculate the expectation value in Eq. (7), we go 
beyond a linear-response model and use the Keldysh- 
Schwinger framework (for a review, see Ref. 24). In such 
a closed time path formalism, the statistical operator is 
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designed to describe electronic transport in a system that 
is out of equilibrium due to an external perturbation 24,25 . 
The operator average takes then the form 




?7=± 


( 8 ) 


The brackets (.. ,)q denote the average with respect to 
the free Luttinger liquid Hamiltonian Hll- 77 = ± is the 
Keldysh index representing the position of the respective 
operator on one of the two time branches of the Keldysh 
contour and Tk denotes the Keldysh time-ordering 
operator 24 . Here, we used f K dt = f^^dt rf. The 
expression can now be expanded up to the desired or¬ 
der in the external perturbation. In first order of the 
impurity Hamiltonian, it can be readily shown that the 
Keldysh time-ordering of the two operators yields their 
retarded correlation function and Eq. ( 8 ) becomes the 
usual Kubo formula 26 . In the following, let us define the 
transmitted current in the system as the difference of the 
free, and the backscattered current, jt T = jo~jbs- In this 
notation, jb s is a positive quantity, reducing the bare cur¬ 
rent due to backscattering off the impurity. Finally, we 
can compute the correction to the conductance (resulting 
from the impurity) by SG = > where we 

defined the correction to the transmitted current opera¬ 
tor as 6j tl = jti - (jtr)o = ~ jbs- 

As an additional quantity, the shot noise at zero fre¬ 
quency indicates the effective charge involved in the 
backscattering process. We calculate the symmetrized 
current-current correlation function 12,25 , which at zero 
frequency reads 

S{u -> 0 ) = ^ f dt 2 {TKSju'it) 2 )). (9) 

v 

The effective backscattered charge, also called Fano fac¬ 
tor, is then given by the Schottky formula 12 e* = 2e \(j b )| 
in the weak backscattering limit. 


C. Effective operators 

The free current operator of the system simply reads 

j 0 = e ^(N + -N_). ( 10 ) 


Since the zero-modes are decoupled form the electron in¬ 
teractions, the bare Fermi velocity vf is implemented in 
Eq. (10). The external bias has the following impact (see 
Appendix B), defining a new, effective free current oper¬ 
ator 


(jo) / 


jo + 


e 2 V 

27T 


In the absence of impurities, we thus find from Eq. (7) the 
averaged transmitted current (j tr )o = ((jo)j)o = 
representing the constant conductance G 0 = of the 
helical edge channel. This contribution is interaction- 
independent, as it should be, which is a direct conse¬ 
quence of the decoupling of the zero-modes from the in¬ 
teractions in this model. It can be seen as a response of 
the system to the bias changing the number of right and 
left moving particles, independently of the system size. 
To account for the influence of the impurity, we next de¬ 
fine the backscattering current. It is given by the rate of 
changing the free current with time due to backscattering 
processes off the impurity 12 . Employing the Heisenberg 
equation of motion, we use 2 ' 


jbs(t) 


i 

2vp 


dx [H R (t),j 0 ]. 


To compute the average backscattering current, we more¬ 
over need to find expressions for the operators (jbs) l 

and > modified by the voltage shift. From now 

on, we specify the form of the impurity. For simplic¬ 
ity, we consider a point-like Dirac-impurity of the form 
a(x) = a ( avF) d(x), where the impurity strength is 
given by the dimensionless parameter 10 a (we take as 
well m[x) = m vfS(x) with dimensionless m in Eq. (5)). 
At this point, the limit L —> 00 can safely be taken, and 
as a consequence, the zero modes N± are neglected. In 
the absence of zero modes, the Klein factors commute 
with all fields and will compensate in the averaged term, 
so we can drop them as well. We derive (see Eq. (B4)) 


.avF f 27ra 


(«*), (,) • „ V L . 

avF ( 27ra N K r 


: (d x 9{0,t) + e 2 i«°-t)+fi) ;+ H.c. 


2t t \ L 


1 

vK 


2 vf J 
')* 


,t) \ . e ieVt 


+ W 

vf 


(Q te ieVt^ +H c _ 


( 11 ) 
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The structure of the backscattering current operator is very similar to the one of the impurity Hamiltonian, but note 
the opposite sign of the conjugate terms 


(jbs) j (t) 


—e 


avF 

7T 



^ 0 ( 0 , 4 ) 



e 2i(0(O,i) + ^t). _ jj c 


= ie 


avp 
27r 



2_: (d t e 2i?i(0,t)N ) :e ieV * 
vK \ J 


J_. e 2icK0,t). (Q te ieVtj _ H c 


( 12 ) 


Since the impurity Hamiltonian is treated perturbatively 
here, we make use of the identity in Eq. (4). Eqs. (11) 
and (12) will then be of the form of a total time deriva¬ 
tive in the cases of vanishing interactions v —> vp and 
K —^ 1, Galilean invariance (vK = vf), or a model with 
fully interacting leads. As we will see below, this cor¬ 
responds each time to a zero backscattering current in 
lowest (second) order perturbation theory of the impu¬ 
rity strength. We observe from Eq. (12) that because of 
(<5jtr)o = — (jbs)o = 0, there will be no correction to the 
current in first order of the impurity strength. The same 
holds for all odd orders as a consequence of the “neutral¬ 
ity rule” for the respective correlation functions 14 . In the 
following, we will evaluate the two lowest order contri¬ 
butions to the backscattering current, that are in second 
and fourth order of the impurity strength. 


III. BACKSCATTERING CURRENT TO 
SECOND ORDER 


A. Bosonic approach 


Having derived the effective operators in the presence 
of the bias, we can now compute the backscattering cur¬ 
rent with the help of Eq. (8). To lowest (second) order 
of the impurity strength, we find (see Appendix C) 


(jbs) 



e ( a \ 2K (AT - v F ) 2 ( v) 2K+i 

2nT{2K) \vJ (AT v F ) 2 1 


(13) 



that in the non-interacting limit the average backscat¬ 
tering current vanishes, as it should be. Importantly, we 
used the fact that the ground state is taken at zero tem¬ 
perature, corresponding to the regime kpT -C eV. The 
Poissonian shot noise can be evaluated analogously from 
Eq. (9). We obtain 


S(u —» 0) — 2e (jbs)- 


The effective backscattered charge in this case is e* = 
2e (? ^ = 1, indicating single-particle backscattering. 
When comparing to the case of the magnetic impurity 


(see Eq. (5)), we obtain in an identical fashion 


(jbs)m — 


27ro 2 T(2K) \v 


/ n \ 2 K 

0 w 


2K— 1 


and the same effective charge. In this case, we observe 
an additional factor of V~ 2 compared to the Rashba im¬ 
purity, due to the lack of the two spatial derivatives in 
the magnetic case. Moreover, while the contribution of 
the Rashba SOC was primarily inelastic, the backscat¬ 
tering off the magnetic impurity does not vanish in the 
non-interacting limit. 

The existence of a finite contribution in second order in 
Eq. (13) relies on the distinction of the two velocities vK 
and vf in our setup. Such a distinction is made at two 
places, in Eqs. (1) and (10), associating vp with the zero- 
modes (modelling non-interacting leads) and vK with the 
bosonic excitations. At this point, we would like to in¬ 
vestigate if such an assumption is in conflict with basic 
symmetries of the system, such as Galilean invariance. 
Indeed, as was originally pointed out by Haldane 16 , this 
symmetry requires vK = vf, e.g. in systems with a 
parabolic energy dispersion. We therefore expect that 
the second order contribution to the current in Eq. (13) 
would vanish if Galilean invariance were to hold. The 
question of Galilean invariance is a well known problem 
of one-dimensional systems. As discussed in Ref. 28, this 
invariance is not automatically accounted for in Luttinger 
liquid theory, since it uses linearization of the energy 
spectrum as a vital requirement for bosonization. When 
dealing with Luttinger liqids in an underlying parabolic 
system, one thus faces the problem of artificially broken 
Galilean invariance 16 . To cure this problem, different ap¬ 
proaches were proposed, either defining new expressions 
for v and K to satisfy vK = Vp 2S , or by a careful treat¬ 
ment of higher harmonics of the density functions 15,29-32 . 
In the scenario described here, however, there is no need 
for such modifications, since Galilean invariance is emer- 
gently broken. Edge states in the gap of a topological 
insulator form a spectrum of two branches with opposite 
slope, that cross at some point in the bulk gap. Except 
for influences of remote bands, those branches are almost 
linear functions of the momentum, such that we find a 
solid ground for bosonization and the application of the 
Luttinger liquid model. The system is quasi relativistic, 
breaking Galilean invariance, and therefore the Luttinger 
liquid relation vK = uf( 1 — g 2 /(27 TVf)) 7^ vp remains in 
general valid. Even though the second order correction 
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to the conductance derived from Eq. (13) is still irrele¬ 
vant in an RG sense, and thus never leads to localization 
in the low energy limit, we find in general a contribu¬ 
tion scaling as SG ~ V 2 in the limit of weak interactions. 
Studies of possible mechanisms to backscatter inelasti- 
cally have revealed various power law corrections, indi¬ 
cating the inelasticity of the process 6,9,11,12 ’ 33 ’ 34 . Among 
such studies, the correction we describe above is one of 
the lowest powers ever found, and is thus expected to be 
of high relevance for transport. 

If Galilean invariance was restored for any reason, or in¬ 
teractions were fine-tuned such that vK = Vf, the SPB 
contribution would vanish, and the lowest order correc¬ 
tion to the conductance would be given by correlated 
two-particle processes, as discussed in the Sec. IV below. 


B. Fermionic approach 


To support our result, we performed an additional 
calculation in fermionic language using the generalized 
Fermi’s golden rule. According to this model, the av¬ 
erage inelastic backscattering current is proportional to 
the transition rate 6,33 \(f\T\i)\ 2 5 eitef , where | i) and |/) 
describe intial and final multi-particle states that, in the 
case of single-particle backscattering, differ by one in the 
number of right and left moving particles, and e/ 
are the energies of the two respective states, and the T- 
matrix operator is given by 26 T = H' + H'GqT. We 
specify the perturbation Hamiltonian H' = Hr + H lnt 
and Gq = e ]; Hq , where Hq is the free kinetic Hamil¬ 
tonian. The average backscattering current is precisely 
defined as 


(jb.) = -Ei%i 2 W// w - ( 14 ) 

vf “ 

Here, M, t f is the transition matrix element, given by 
Mjf = (f\T\i). We focus on single-particle contributions, 
where M,/ is of order a <72■ Therefore, we go to second 
order in T = H'GqH' and consider the terms 

M if = (f\H'G 0 H'\i) = (f\H R G 0 H int + H int G 0 H R \i). 

(15) 


We next derive the form of the effective operator 
T in momentum space, where we define for the 
Fourier components 35 T^x) = Efc e +lkx ^ ± k . 
The free Hamiltonian can thus be written as Hq = 
Efc £+(fe)^+ )fc ^+,fc + £-(fc)^L ifc ^-,fc with single-particle 
energies £±(k) = for the right and left movers. 

The impurity Hamiltonian in momentum space becomes 


Hr 


J dx a(x) (E’I'V) - T f + (x) +H.c. 


E r„ 6 (^ ifca $_ lfcl - tf_ iki 9 + , ka ). (16) 

k a ,kb 


with a(x) = aLvF S(x) and T a b = iavp(k a + kb). The 
interaction Hamiltonian on the other hand reads 

Hmt = ~9i J dx E (x)E(x)'I' + (x)'I , _(x) 

= Z f L E ( 17 ) 

k r ,k s ,q 

In this notation, a is a dimensionless parameter, while g 2 
has the dimension of a velocity. Now, we can compute the 
two parts of T in Eq. (15). Using the anticommutation 
relation = {4b, Tj} = 6 itj = S kitkj Sr, itVj , 

we find after some algebra 


T\i)=P e 

k a ,k s ,q 




_,_ fea+9 ^ + , fes+ g* + ,fc s *+,fc, 






(18) 


with /3 = ~ 8 ” g2 . Here, we made use of energy conserva¬ 
tion, requiring — (k a — q) + (k s + q) — k s — kb = 0, meaning 
that by setting kb = 2 q—k a , the sum over one momentum 
can be dropped 36 . As expected, the energy denominator 
cancelled with the momenta in r a ^, regularizing the ex¬ 
pression. 

Let us briefly interpret Eq. (18). Thinking of 

Efc s ’^ , ±,fc s ±<z’^±>'b> = P±.±g 85 a P ar ticle density, the 
effective operator T exhibits the same structure as 
the bosonized Rashba Hamiltonian in the first line of 
Eq. (11). It describes a backscattering process coupled 
to a bosonic excitation at one of the branches and, im¬ 
portantly, implies finite interactions. Without interac¬ 
tions, we can readily check that the transition matrix 
elements (f\H R \i) and (f\H R GoH R \i) are always zero be¬ 
cause of energy conservation. This result embodies the 
non-interacting limit and was shown already in Ref. 10. 
The components with q = 0 in Eq. (18) can be identified 
with the term (A+ - iV_) e 2 ^(°) ~ Ve 2 ^ in Eq. (11), 
and thus manifest the zero-mode part of the T-matrix. In 
the bosonic calculation, we assumed that the zero-modes, 
representing the external contacts, were decoupled from 
the electron interactions. To adopt the same model here, 
we have to impose an additional constraint on the effec¬ 
tive operator in Eq. (18) in such a way, that the inter¬ 
acting, effective operator does not alter the zero-modes. 
This can be achieved by excluding the q = 0 compo¬ 
nents explicitly from Eq. (18). Indeed we find from an 
explicit calculation (see Appendix D), that Mif = 0 if 
all values of q are taken into account in Eq. (18). This 
case corresponds to the situation when interactions do 
affect the contacts, and the finding is thus in agreement 
with Eq. (13). On the other hand, if we exclude the 
components with q = 0, we get Mif ^ 0 and a finite con¬ 
tribution in second order of the impurity strength. This 
case reflects the setup we consider in this article, with 
the leads decoupled from the electron interactions. 
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Defining a system with a finite, effective bandwidth, we 
numerically evaluate Eq. (14). Hereby, we fix |i) as an 
initial state of the system in the presence of an exter¬ 
nal bias, and sum over all final states |/). In both the 
initial and the final state, we account for a Fermi sea of 
finite depth, and finally estimate the dependency of jbs 
on the voltage bias (see Appendix D). We conclude, that 
the backscattering current can be approximately given 
by the expression 

(jbs) ~ ea 2 g | 4 ^ (eH) 3 , (19) 

which is in qualitative agreement with the result in 
Eq. (13), in the limit of weak interactions. 


IV. BACKSCATTERING CURRENT TO 
FOURTH ORDER 


The analytical calculation of the average backscatter¬ 
ing current in fourth order is a challenging task, since 
in the Keldysh scheme, multi-dimensional coupled inte¬ 
grals are encountered. In principle, it will turn out to 
be very tedious if not even impossible to account for all 
the contributions arising. Such contributions can quite 
generally be processes involving two particles that are 
coupled by interactions. For simplicity, we consider here 
only the two extreme limits of the problem. Those are 
the cases of strongly coupled two-particle backscattering 
and the backscattering of two decoupled particles. Tech¬ 
nically, these limits are taken by performing contractions 
on the time variables. The remaining integrals can then 
be solved to reveal the respective contributions to the 
backscattering current (see Appendix E). 

We first focus on the strongly coupled TPB limit. It is of 
particular importance, since TPB at strong interactions 
was shown to be a potentially relevant perturbation in 
the RG sense, leading to localization and a breakdown of 
the conductance in the low energy limit 9,10,13 . qq ie TPB 
contractions are only meaningful if the processes selected 
in the procedure represent the most important contribu¬ 
tions of all processes. In particular, this requires that the 
time variables we consider as small in the contraction pro¬ 
cess appear indeed with sufficiently fast decaying power 
laws. For the Rashba impurity, this is the case only in 
the regime K < 1/2. For such strong interactions, we 
can infer that the TPB processes are of leading power in 
the bias and we find (Appendix E) 


(jbs) 


7t2 8A+5 ea A 
K 2 T(8K) 4(2t r) 4 


(°r 2 


(v)V> 


8K-1 


( 20 ) 


if K <1/2. Likewise, the noise can then be derived as 
S(u -> 0) = 4e(j bs )- 


We hence obtain an effective charge e* = 2. This result, 
in addition to the expected scaling of V SK ~ l , indeed con¬ 
firms a pure TPB process 9 ’ 10 " 12 ’ 37 . 


Studying the same contraction procedure for the mag¬ 
netic impurity, we see that these TPB contractions are 
never justified, because of the lack of derivatives in this 
case. Therefore, terms as found in Eq. (20) are never the 
leading contributions. The TPB current given above sur¬ 
vives in the case of restored Galilean invariance, vK = 
Vp, in contrast to the SPB processes of Eq. (13). In an ex¬ 
periment at K < 1/2, we should find e* = 2 from TPB in 
the presence of Galilean invariance, while in the absence 
of this symmetry, the second order SPB contribution was 
dominant, leading to e* = 1. Consequently, besides the 
distinct power scaling, we can identify the Fano factor as 
a direct evidence for the presence of Galilean invariance, 
which should therefore be experimentally measurable. 
Next, we consider the situation of weak interaction 
strengths K > 1/2. In this regime, the TPB contrac¬ 
tions used before are not well justified, since they do not 
capture the essential subset of processes any more. How¬ 
ever, we can still observe the change of the scaling of the 
strongly coupled TPB terms. When going from K < 1/2 
to K > 1/2, we infer from the structure of the encoun¬ 
tered integrals that the scaling of such terms is changed 
from V 8K ~ 1 to V 4K+1 . We have to keep in mind though, 
that the leading contributions for K > 1/2 are not rep¬ 
resented by pure TPB processes any more, but should be 
more of the form of two weakly coupled backscattering 
events. Along this line, we can confirm the existence of a 
crossover of scales at K = 1/2, for correlated two-particle 
processes, as discovered in Ref. 10. 

Finally, to complement our fourth-order analysis, we con¬ 
sider the opposite limit of two processes decoupling into 
separate single-particle events. Here, we perform the con¬ 
traction of time variables in a different way to sort out 
the terms of interest (see Appendix E) . We then find 


(jbs) 


(vK ~ vf) 4 a 4 e /«\ 4A+1 ( y) 4K+2 

(vK) 4 167r 2 r(2K) 2 \vJ 1 ’ 

( 21 ) 


S(uj 0) ~ 4e(j bs ). 


In this case, the contribution has a similar structure as 
the backscattering off a magnetic impurity, where we can 
derive 


(jbs) r 


me 


(“f 3 (?)'<-> 


167T 2 r(2K) 2 \V 


4K—2 


( 22 ) 


and the same Fano factor. Since two times a single charge 
is transfered by the two decoupled SPB events, we find an 
effective charge of two, meaning that we can not distin¬ 
guish double SPB and TPB from the noise only, but from 
the scaling with the bias. This kind of SPB contraction is 
always justifiable, which is why the current contribution 
in Eq. (21) holds for all 0 < K < 1, however, the TPB 
terms will still dominate for strong Coulomb interactions. 
Consistently, the SPB scaling of U 4A+2 exhibits exactly 
one more power of V than the above mentioned corre¬ 
lated two-particle processes for K > 1/2. The missing 






power in the latter case reflects the existence of one ad¬ 
ditional time integral, creating a weak link between the 
two events. 

With Eqs. (20) and (21), we have given the two inter¬ 
esting limits in this order of a, being strongly coupled 
two-particle backscattering and decoupled double single¬ 
particle backscattering, respectively. The intermediate 
regime can not be captured exhaustively in our present 
analysis, though a transition between the two limits is 
indicated by the power scalings found from the structure 
of the integrals. 

V. CONCLUSION 

In this article, we have studied the influence of a local 
Rashba spin-orbit scatterer on the edge transport of a 
helical quantum spin Hall system. In contrast to previ¬ 
ous analyses, we have applied a non-equilibrium approach 
with explicit implementation of the contacts attached to 
the sample. Indeed we have found that the modeling of 
the leads is of qualitative importance for the resulting 
backscattering current in lowest order. Because an emer¬ 
gent Lorentz invariance is characteristic of the QSH edge 


channels, we have discovered a single-particle backscat¬ 
tering allowed at second order of the impurity strength. 
This contribution differs from the backscattering off a 
TRS-breaking impurity by its purely inelastic character 
and by the power law scaling. 

Subsequently, we have extended our calculation up to 
fourth order in the impurity strength, to verify the exis¬ 
tence of correlated TPB processes. As expected we have 
identified the corresponding power law contributions of 
the correction to the conductance, as well as the crossover 
of scalings at K = 1/2, in agreement with Ref. 10. The 
analysis of the shot noise has yielded the backscattered 
charges in the limits of pure SPB or TPB. The resulting 
Fano factor of one or two, respectively, can be seen as 
a direct evidence for the absence or presence of Galilean 
invariance. 
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zwerk Bayern 11 (ENB graduate school on “Topological in¬ 
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Appendix A: Form of the interaction Hamiltonian 

Let us comment here on the specific form of electron-electron interactions in a helical liquid. In a usual, spin¬ 
ful Luttinger liquid, Coulomb interactions are generally described by two-body interactions of the form = 

J dx f dy (cc) + (z/) C7 (a; — y) (y) HI (a;) 16, 17 > 26 > 38 , with operators arranged in normal order. H/^a;) and T(x) are 
fermionic creation and annihilation operators for an electron at position x, that will be specified with additional in¬ 
dices for momentum and spin. In general, at the QSHE edge, we expect only short range ineractions due to screening. 
Assuming contact interactions, U(x—y) = Uq S(x — y) , we can classify three types of interactions H int = Hi + H 2 + H 4 , 
with 

H l =U 0 ^2 [ dx 

ct,<t'=± 1/2"' 

H 2 = Uo ^2 [ dx ^+,Ax)^-y( x )^-,a'(.x)^ + ^(x) 1 

<7,<T'=± 1/2"' 

Ha = U 0 J2 [ dx y^A X )^ + ,*>( X )' ^/ + ^( X ) i^f +,°( X )■ 

Here, we have specified the fermionic field operators such that a (x) creates a right (+) or left (—) moving particle 
with spin a = ±1/2. In the helical liquid, because of spin-momentum locking, one index is redundant and can be 
dropped again later. In terms of the usual g-ology 14,28,30 , we could replace, in each of the Eqs. (Al) to (A3), Uq by 
prefactors Si,||/_l = S 2 .||/_l = 9i,\\/x = where the index || indicates parallel spins a = a 1 and ± represents a = —a'. 
Here, the Pauli principle is manifested in forbidding terms <? 4 ,||, since in Eq. (Al) ip +t<j (x)ip+ t( j(x) = 0. Moreover, 
we see that after fermionic anticommutation gi j| is of the same form as g 2 ,|| but with a minus sign, and therefore 
cancel each other in case of the contact interactions we assumed here. Because of this, all the parallel interactions are 
absent. 

In the helical liquid, we have a special situation because of spin-momentum-locking. This tells us, that processes 
<74 ,-Lj 3 i,_l (and gi.\\, g 2 ,\\) are impossible due to the helical character of the edge states, therefore, the only remaining 
interaction term is g 2 ,±, which is 

H 2 = u 0 f dx ^+,a( X )^-.-a( X )^--<r{ x )^+,cr{ x )- 

< 7 = 1-4 J 


(Al) 

(A2) 

(A3) 
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Since spin and momentum are locked, one of the two indices can safely be dropped. Importantly, the helical liquid is 
only formally spinless, but as we pointed out, there is an essential difference between the helical and the spinless case, 
when it comes to Coulomb interactions. In a spinless picture, we would expect no contribution at all in case of contact 
interactions, since 34 cancels due to the Pauli principle and as well 52—31 = O 23,30 . From these considerations, for 
contact interaction, we have in the helical liquid 

Hint =52 J dx ^(a:)^(a:)^_(x)^ + (a:). 


Appendix B: Derivation of the shift of bosonic fields 


Starting from Eq. ( 6 ), we rederive here the shift of fields as presented in Eq. (7), following schematically Ref. 19. 
To do so, we introduce an unitary operator U, defined as 

U = e^n/W-). (Bl) 

a is a constant factor not yet specified. Here, the Klein factors are written in exponential form, fj. = e l H and 
F± = e~ l H. Those Klein factors commute with the bosonic creation and annihilation operators, though not with 
the particle number operators. The commutation relations read [N v ,/„>] = —i6 v „>, and [f v ,fri'] = [N V ,N V >] = 0 with 
77 = ± i8 > 19 . We now see that U is able to establish the shift 

= e -0(H LL -H v ) e c = e -p H’ Ve c, 

when a = eL/(AvF^) is chosen. e c = e i ^ v “ vv f/ l is just a constant created by commutating U past ef iHv . Applying 
U t from the left hand side is equivalent to a voltage-dependent shift on the particle number operators N± —> N± ± . 

With this, we write now Eq. ( 6 ) as 


(O) = -Tr {Ue~ pHhL U^ e iHt Oe~ iHt ) 

Zll 

Tr(e~^ HLL i{HLL+Hy-\-Hii)t j 


Zll 


(B2) 


with Zll = Tr(e~^ HLL ) and H = Hll + Hr. Moreover, we define for a general operator A = ( WAU ). Above, 
we used the cyclic permutation of the trace and inserted factors of UU' = 1. Since the system is mainly governed 
by the exponents linear in time, we approximately neglected all the terms of the power of t 2 , leading to U'e lHt U ~ 
e i(H LL +H v +H R )t e C/p _ Next, we go to an interaction picture 19 , using that for any operator Ai(t) = e lHvt A e~ lHvt , 
we write 


e -it(H LL +H R +H v ) _ ex p 


—1 


dt H LL (t') + (HR). f (t / ) 


= e~ itHv S{t). 


Eq. (B2) then becomes 


{(D) = -J— Tr {e-P HLL S\t)e itHv 6e- itHv S{i)) = -J— Tr {e~ pHLL S^(t) (o) (B3) 

Zll Zll k > 1 

and we arrive at Eq. (7) of the main part. In brief, the result in Eq. (B3) means, that introducing a finite voltage of 
the above form is equivalent to transforming both the perturbative impurity Hamiltonian Hr —> (Hr'J (t) and the 

observable under consideration O —> (t). 

The shift of the bosonic fields now depends on the form of the operators appearing in the system. In our case we 
find -A <j>{x) + and d x 6{x) -A d x 0{x) + ^ . (cf. Ref. 12). Performing the shifts on the Rashba impurity in 

Eq. (3), we find with the shortcut 7 = ^ (- jf L ) K 

(Hr^ {t) = F\f _7 J dx 1 : (d x 6{x) + N+ - N-) + e 2 d<K*,t)-^(N + +N_)+^t). + H c ( B4 ) 
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Appendix C: Backscattering current to second order 


Here, we explain in detail, how to exploit normal-ordering to derive the average backscattering current in Eq. (13). 
First of all, to keep the notation compact, we introduce the scalar v~ l (dt) with = 

if dt acts on the bosonic part and v~ 1 (d t )dte leVt = (vF)~ 1 dte leVt , if dt acts on the voltage 
part. With this, we can write Eqs. (11) and (12) in a compact form 


(^«) J (t) = ^ (v-Hdt) dt {-.e 2i ^:e ieVt ) + H.c.) , 

(jb »)j ( t ) =ie^~ (v-\d t ) d t (-.e 2 ^:e leVt ) - H.c.) . 


(Cl) 

(C2) 


To compute (jbs), we start from Eq. ( 8 ), 


Obs it)) = O’bs)/ (*)(“*) J R dt2 (^) 


= PlOL 


2 \ c 

J dt 2 9(t - t 2 )u~ 1 (d t )v'" l (d t2 ) 


(-.d t e 2i ^°'d+ ieVt : - H.c.) , (:d t2 e 2i ^° M)+ieVt2 : + H.c.) 

= Pi a 2 J dt20{t-t 2 )v~ 1 (dt)v~ 1 {dt 2 ) 

x (d t d t2 (h[t.-t 2 )e ieV{t ~ t2) ^ (:e 2imt) e- 2i,p ( 0 ' t2) :^-dtdt 2 (h(t - t 2 )e~ ieV(t ~ t2) J ^ :e -*Ho,t) e 2i4{o,t a ).^ 

+ d t d t2 (h{t 2 - t)e~ ieV(t - t2 ^ ( :e 2 i<Ko,t 2 ) e - 2 imt).^ _ g t g^ ( h ^ 2 _ t yeV(t-t 2 )^ ^ ;e -M*(o,t a ) e 2 i*(o,t).^ 

J dr 0 (r)v~ 2 (d T )d 2 ^h(r) sin(eW) + h(—r) sin(—ekr)j . 


= (-2i)pia' 


/ 27ra 


(C3) 


Here, we used p\ = e and the normal-ordering relations :e 2l ^d) ::e 2 i<t>(t 2 ). _ ;e 2 uf,(t) .. e 2 uj>(t 2 ). _ ^ _ 

t^-.e 21 ^) e -‘ 2l H t 2 )■ w ith h{t\ — t 2 ) = (2ir/L(iv(t\ — t 2 ) + a))~ 2Kl8 . Importantly, time derivatives acting on the above 
operators will create normal-ordered operators in non-exponential form, that will again vanish under averaging, since 
per definition (:A:)o = 0, but (:e A :)o = 1 for any bosonic operator A 39 . This is, however, only true, if the expectation 
value is taken with respect to the zero temperature ground state. Indeed, the expression simplifies a lot, exploiting that 
(:e ±2l ‘^( tl )e :F2l ^( t2 ):)o = 1. After the substitution r = t — t 2 we find the expression in the last line of Eq. (C3). Using 
the fact, that the functions appearing in the integrand only depend on the time difference r, we write d t2 = —d t = —d r 
to simplify 


(jbs) = (-2*)pia 2 J drv 2 (<9 r )<9 2 ^/i(r) sin(eUr) 

= 2 ipia 2 yi v (eV) 2 j drh(r)MeVr). 


(C4) 


In the last step we applied the definition of the v 1 (d T ) and performed an integration by parts. Evaluating the integral 
we arrive at Eq. (13). 


Appendix D: Fermionic approach 

For practical reasons, we first rewrite T in Eq. (18) with positive momenta. Considering a system with only positive 
energies, we can transform k —> —k for all left-moving states, for convenience. The states are then labeled by their 
discretized momenta fcj = ^n,, where the integers n, = 1 , 2 ,... 2 tifs are limited by the parameter n^s, and we find 
now for both species £(£',;) = vpki = Vf j/ n i- The energy of the multi-particle state is the sum of the energies of all 
contained single-particle states ’E-t*,.. We now define the initial state | i) as follows: A Fermi sea of finite depth is 
represented by ufs occupied states on both branches. Next, we apply a voltage bias in such a way that the ±-branches 
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Log(p) 



Figure 2. Log-Log Plot of the numerical factor p and the voltage integer Nv {Nv = Tj 4 ^)■ The continuous line represents 
p = (21Vy ) 3 for comparison. 


in the presence of the bias are filled up to rii < hfs ± Ay = Hfs ± with Ay < Here, by introducing the 

integer Ay, we made the choice of applying the bias symmetrically on both branches. We thus have 


I*) 


riFS+Nv 

n 


3 =1 



riFS-Ny 

n 


3 = 1 



| 0 ), 


where |0) is the vacuum state. Hereby, we defined a (arbitrary) reference order of the right and left moving single 
particle operators. This order has to be kept throughout to keep track of all signs. The two parameters to be set 
in this model are ufs and V, where ufs can be considered as a symmetric momentum cutoff. For a general set of 
nps and V, we note the following. First, as explained in the main text, the matrix element M; t f = 0 if all values of 
q are taken into account. The cancellation mechanism is, that one term with nonzero q will always be compensated 
by another one with q = 0 in any of the terms given in Eq. (18). If we exclude the components q = 0 from Eq. (18) 
though, we get M-if ^ 0. Second, for a positive bias, the processes backscattering a left mover into an additional right 
mover, that are represented by the {+ o —} in Eq. (18), are always zero (and the other way round). Third, the 
effective cutoff (or depth of the Fermi sea taken into account) nps does not matter. Because of various cancellations, 
only processes close to the Fermi surface contribute to the backscattering current. 

Adopting the same scheme as in the bosonic part, with the zero-modes decoupled from the electron interactions, we 
exclude the elastic q = 0 components from Eq. (18), and subsequently employ Eq. (14) to derive the backscattering 
current. Scattering processes leading to the same final configuration will be added with their respective signs, that 
arise by sorting the single particle operators of the many particle state in the reference order. Adding all contributions 
according to Eq. (14) yields a current 


Obs) = — 

Vf 


-iag 2 


2 2 
ea J g$ 

V = - yP- 

Vf-L 


Here, we have introduced the numerical factor p, representing the sum of all diagrams. In general, p will increase 
with increasing bias. Realizing that the result does not change with the momentum cutoff, the only parameter of the 
system is just V. 

To estimate the dependency of p on the voltage V, we can fit a power law to the numerical data, where we would 
expect from the bosonic calculation (see Eq. (13)) jb s ~ V 2K+1 ~ V 3 for weak interactions. Indeed we see that 
a scaling of p ~ (2Ay) 3 = (2^|^) 3 is a reasonable approximation, particularly for large voltages (see Fig. (2)). 
Deviations from this scaling for small values of V are probably due to finite size effects of the calculation. We finally 
arrive at Eq. (19) of the main text 


L 2 

v%{2f) 3 


(eV) 3 . 


O'bs) ~ ect <?2 


(Dl) 
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Appendix E: Backscattering current to fourth order 


In this section, we give details of the calculation in fourth order of the impurity strength. This means, we expand 
Eq. (8) to third order in the impurity Hamiltonian 


7 )=± 


(jbsfa)) = ^ E \^ K (• ?bs )/ ^6 J dtzdtidts i (t 3 ) (Hr) T fa) ^HrJ^ fa) 

/»00 

/ dt 3 dt\dt 3 r( rf' rf" 


J / o 


_ P2_ ( 27 ]a\ iK nA 


«* E 

rj,rj',i}",rj"'=± J 00 


12 \ L ) 

x( t k Y, ( C S!wW" ^ 3 ’ 4 ’ 5) - -»■ "^{ 2 , 3 , 4 , 5 }) ) ■ 

' a 3 ,a 4 ,a 5 =± / 0 


(El) 


With p 2 = / 2 ^)4 ■ Here we have defined a shorthand notation 


cffiXu 2,3,4,5) = v-H$>~ 1 ($ 3 )«'~ 1 ($>~ 1 ($.) 


x d t2 dt 3 d t4 dt 5 


(h(2 , 3)h(2,4 )/i(2,5)/i(3,4)/i(3,5)/i(4,5)) e ie WA2t2+A 3 t 3 +A 4 t 4 +A 5 i 5 ) 


x . e iA 2 20’>(2) e iA 3 20'' (3) e iA 4 2^ (4) e iA 5 20” (5). 


(E2) 


A* = ±1 just represents a sign here and hfa,tj) = (2n/L(iv(ti — tj) + a)) 2KXiXj with i,j £ {2,3,4,5} 18 . The 
neutrality rule 14 holds overall, such Xu =2 = 0- From the same arguments used in the previous section, we argue, 
that the derivations act only on the terms in front of the operators in normal ordering signs at zero temperature. 

We can see from Eq. (E2), that (2,3,4. 5) stays invariant under a simultaneous change of two times ti and 


tj, A i and A j (and the Keldysh index r/i and ijj). Importantly, this is only true if the order of operators e lXi ^ l ' ) can 
be interchanged, which is possible before applying Tk , but not after. From this fact we infer, that we can change for 


instance C+ J) ,+ „ ?7 , / 


(2, 3,4, 5) and C+~r+^ 


(2, 3,4,5) to C++, v „ v „,( 2, 3,4,5) in Eq. (El), since the time variables 


t 3 , ^ 4 , t 3 are integrated over and can just be renamed (this brings a factor of 3 in front). We get 


{jb s fa)) = Y a *(~T-~) E [ dt 3 dt i dt 5 r] , r]"r]' 
4 \ h J ViV , iV ,, tV ,„=±J -°o 


x ( T k { C tvWW- (2,3,4,5) - C--+t v ,„ (2,3,4,5)) 


(E3) 


The next step will be explicit time-ordering. Because multi-dimensional, coupled integrals are involved, we 
were not able to find an analytic result for the average backscattering current. However, to proceed, two different 
approximations are studied, in the form of contractions in time. Physically, we assume, that among the four operators 
with four different times occuring in the integrand, some operators couple strongly to each other, while others only 
couple weakly. This procedure is mathematically justified under certain conditions, that we specify below. 


Time-ordering We continue from Eq. (E3) and split the full expression into two parts where the one depends 
explicitly on the Keldysh indizes, and the other does not. 


/ c\ \ 4 K /»oo 

(jbsfa)) = J dt 3 dt A dt 5 C%’ v ,v ’ v fa,t 3 ,u,t 5 ) (Cffa,t 3 ,U,t 5 ) - Cffa,t 3 ,t A ,t 5 )) . (E4) 
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We undid the full normal-ordering to apply Wick’s theorem in the following. Now, the terms are 


C™ , ' ,1 "’ r ''"(t 2 ,t 3 ,t 4 , i 5 ) = £ vW'v^lTK^Wr.e^'W:: 

77 , 77 ' , 77 " ,rj'"=± ' 


e -i2r (4) ;;e -<2^ (5). 


E v'v"v'"(T K : 


-t2^(2) ::e -i20*» (3).. p +t2^ (4).. p +i2^» (5). 


'::e 


::e 1 


77 , 77 ' , 77 " ,77 ,// =dz 

Ci(t 2 ,t 3 ,t 4 ,t 5 ) = v~ l {d t2 )v~ 1 {d t:i )v~ 1 {dt i )v~ 1 {dt b ) 


( (d t2 dt 3 d u d t5 ') 

(h( 2,3)/i(2, 4)h(2, 5)h(3, 4)h(3, 5)/i(4,5)) e ±ieV(t a +t 3 -u-t s ) 

\ 

V 

(h( 2,3)/i(2,4)/i(2,5)/i(3,4)/i(3,5)/i(4,5)) 

) 


/ // /// 

Cq' v ' v ' Jl (t' 2 -t 3 , £ 4 , < 5 ) is the same for both configurations, because of the Debye-Waller formula 14 . Since the average 
acts only on operators, we can pull out the factor including the time derivatives, keeping in mind, that this factor 
depends on time ordering, indicated by the label | or d- However, it is still a scalar, which is why we can use Wick’s 
theorem to simplify the product Co of the exponential operators. 

The next step is to perform time-ordering explicitly, so to evaluate the sum over Keldysh indizes. Hereby, we use 
Wick’s theorem, which for exponential operators is nothing but the Debye-Waller relation 


Cr ,V " (t2,t 3 ,U,t 5 ) = V W" II (T K :e i2X *+*&h:e i3X *+ ,,i W: 

v,v'iV "!»?" , =± i<j,(i,i)e(2,3,4,5) ' ^ 0 

= E W ?5 Cq 2V3 ( 2 , 3)Cq 2V4 ( 2 , A)Cq 2115 ( 2 , 5)Cq 37?4 (3, 4)Cq 3V5 (3, 5)Cq 4T?5 (4,5). 


ri2,ri3,V4.,V5=± 


(E5) 


For now, rji, rjj £ ( r/ 2 ,7?3,774,775) have to be identified with the (77,77', 77", 77'"). We defined 

Cr J (x,y) = (T K :e l2Xi ^ r ' z ^::e l2Xj ^'' J ^:} o- At this point, we have to explicitly perform the time-ordering. To do so, 
we specify: 

Co + ( x , y ) = y)°( x -y) + Hv, x )Q(y - x ), 

c£~(x,y) = h(y,x), 

Cq + ( x i y) = h{x,y), 

c o —( x >y) = H x ,y)8(y ~ x ) + h(y,x)6(x - y). 

With this, we perform the sum in Eq. (E 5 ), eliminating along the way mute combinations with unphysical, cyclic 
time-ordering. 

TPB Contractions In this subsection, we contract times with a focus on the possibility of TPB. Those are processes, 
where two particles are backscattered from left-movers to right-movers in a very short sequence, and the reversed 
process takes places at another time. In the contraction process, we use the abbreviations y = t 2 — t 3 , y' = £4 — t 3 , 
Y = {t 2 + t 3 )/2, Y' = (t± + t§)/2. Then, we replace the integration variables t 2 ,t 3 ,ti by y,y',Y’. The contraction 
implies, that y, y' « Y,Y', so we neglect the y, y' whereever safely possible, treating the cutoff a with care. 

We now Taylor expand the expression around the small times y, y' ~ 0 , and find terms of different powers of the times 
y,y' and (Y — Y'). In agreement with the spirit of the contraction, we focus on leading powers, that are 21\ — 2 for 
the small times y, y', and —8 1\ for the time (Y — Y'). From these powers we can infer, that the TPB contraction 
is only well-justified for K < 1/2, since only then the functions containing the small time variables y, y' decay fast 
enough. On the other hand, for K > 1/2, the TPB contributions selected by this contraction procedure, still exist, 
but will be dominated by other terms where the y, y' are not small time distances any more. 

Proceeding as described above, we find in lowest order the full term 


Cq ' v ,i? ,v (t 2 , t 3 , £4, £ 5 ) (Ci(t 2 , t 3 , t±, t 3 ) — Oi (t 2 , t 3 , <4, t 3 )) 

= pvp [ T ) (- 16 )(! - 2 KfK 2 v\ Sl S2 )( Si (L 2 L 3 ) + s 2 {L\ - La)). (E6) 
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Eq. (E6) can now be applied to Eq. (E4), after evaluating the three different types of integrals: 

Li/a ~ J dY(TivY + a)~ 8K 0(Y) ( e ~ 2ieVY - e +2ieVY ) = - i2 8K e ±4iK7T cos(4ATr)F(l - 8 K)(eV/v) aK - 1 v~ 1 3 , (E7) 


La 


dY((vYy)-* K 0(Y) (e 


— 2 ieVY + 2 ieVY\ 


— e 


) = —i2 8K cos(4ATr)r(l - 8K)(eV/v) 


8 K ~\r\ 


Hereby, we substituted for simplicity Y — Y' —> Y . The short time integrals read 


,2/f-l 


0 ±.iTvK 


s l/2 


(vx) 


2K-1 


dy{±ivy + = =F + ([ _ 2R)y 


(E8) 


(E9) 


In Eq. (E9), the integration range was limited by a general threshold x, which is introduced, since in the contraction 
process, the times y, y' are assumed to be small. Moreover, we approximated x a. Eq. (E9) indicates, that the 
times we consider as small in the contraction process are described by a sufficiently fast decaying power law only in 
the regime K < 1/2, so in this case the TPB contraction is justified. 

For such strong interactions, the TPB processes are of leading power in the bias and we can find from Eqs. (E4) and 
(E6) 


O'bs) 


tt2 8K+5 ecu 4 5 6 

K 2 T(8K) 4(2tt) 4 Vu 




8K-1 


if K < 1/2. 


(E10) 


As a nice verification, the very same result can as well be achieved by integrating Eq. (E4) by parts. From all the 
possible terms arising, we can then source out the only one of the power V 8K ~ 1 , that is exactly Eq. (E10). 


At the point K = 1/2, the contraction process stipulates, that the average current vanishes (see Eq. (E9)). 
For the case K > 1/2, the TPB contractions are not well justified any more. Assuming that we still follow the same 
subset of terms, the rr-dependent term in Eq. (E9) will be of importance. We argue, that the time integration is then 
restricted by the next greater energy scale, being x —> 1/eV. By this virtue, we derive the scaling of the backscatter- 
ing current of ~ V 4K+1 . This means, that we expect a crossover of scales from y 8A_1 to y 4K + 1 at the point K = 1/2. 


SPB Contractions Alternatively to the contractions described above, we can also aim for SPB contributions. 
Hereby, we contract other time-variables together, that come with opposite signs in the exponents. For terms C-, , 
as they occur in Eq. (E4), we define for instance y = t 2 — t< i, y' = t- ti — t 5 , Y = {t -2 + t±)/2, Y' = (f 3 + i 5 )/2, and 
approximate y,y' « Y, Y'. This corresponds to the case where two SPB events are separated in time by a large 
distance, and thus decouple. It turns out, that all the possibilities to contract in this way have to be added up to 
come to a meaningful result. We find 


(jbs) 


('vK - VF) 4 a 4 e (a\* K+ i (pV \4,k+ 2 
{vK) A 167 t 2 T(2A") 2 \vJ 1 ’ 


(Ell) 


Interestingly, this scaling of a 4 y 4A '+ 2 in Eq. (Ell) is achieved in the following way: After the contractions, terms 
proportional to (V— Y') cancel altogether, while from the time-ordering we are left with some Heaviside-functions to be 
evaluated f dY9(Y — Y')Q{Y' — Y) ~ 1/4-. It can be seen, that with this approximation, the correlations between the 
backscattering processes are cut, leaving behind two completely independent SPB events. Those, however, always have 
the same “center of mass”-times Y = Y'. Because of the mentioned decoupling, one time-integral can not contribute 
a factor of 1/P, increasing the scaling by one. Another option would be to contract three times together, which is 
expected to yield SPB contributions of the scaling a 4 V 2K+1 , adding corrections to the second order expression. 
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